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ABSTRACT 

We present a study of the properties of the star-cluster systems around 
pseudo-bulges of late-type spiral galaxies using a sample of 11 galaxies with 
distances from 17 to 37 Mpc. Star clusters are identified from multiband HST 
ACS and WFPC2 imaging data by combining detections in 3 bands (F435W and 
F814W with ACS and F606W with WFPC2). The photometric data are then 
compared to population synthesis models to infer the masses and ages of the star 
clusters. Photometric errors and completeness are estimated by means of artifi- 
cial source Monte Carlo simulations. Dust extinction is estimated by considering 
F160W NICMOS observations of the central regions of the galaxies, augmenting 
our wavelength coverage. In all galaxies we identify star clusters with a wide 
range of ages, from young (age < 8 Myr) blue clusters, with typical mass of 10^ 
M0 to older (age > 100-250 Myr), more massive, red clusters. Some of the latter 
might likely evolve into objects similar to the Milky Way's globular clusters. We 
compute the specific frequencies for the older clusters with respect to the galaxy 
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and bulge luminosities. Specific frequencies relative to the galaxy light appear 
consistent with the globular cluster specific frequencies of early-type spirals. We 
compare the specific frequencies relative to the bulge light with the globular clus- 
ter specific frequencies of dwarf galaxies, which have a surface-brightness profile 
that is similar to that of the pseudo-bulges in our sample. The specific frequen- 
cies we derive for our sample galaxies are higher than those of the dwarf galaxies, 
supporting an evolutionary scenario in which some of the dwarf galaxies might be 
the remnants of harassed late-type spiral galaxies which hosted a pseudo-bulge. 

Subject headings: galaxies: star clusters; galaxies: evolution; galaxies: bulges; 
galaxies: spiral 



Introduction 



Over the last two decades, observational evidence has been accumulating that not all 
spiral galaxies possess a bulge resembling a small elliptical galaxy with a de Vaucouleurs R^/^ 
light profile. Bulges with a surface-brightness pr ofile that does no t follow the de Vaucouleurs 
law also show disk-like, cold stellar kinematics (jKormendylll993l) and light prof i les th at are 
well- fitted by an exponential or an intermediate Sersic law (lAndredakis et al.lll995l ). For 
these reasons these bulges are known as pseudo-bulges or exponential bulges. HST imaging 
of pseudo-bulges ha s also revealed tha t they often contain central compact sources, most 
likely star clusters (jCaroUo et al.l Il998l ) and that they have shallowe r nuclear slopes than 
traditional bulges with the same magnitude (ICaroUo &: Stiavellilll998l ). It is also intriguing 
that the nuclear cusp slopes of the exponential bulges are similar to those of dwarf elliptical 
galaxies with the same radii and luminosities. This similarity suggests that there might be 
an evolutionary link between pseudo-bulges and dwarf ellipti cals, as if the presenc e of a disk 
did not affect the nuclear properties of pseudo-bulge spirals (IStiavelli et al.ll200ll ). 



In order to study in detail the properties of pseudo-bulges and investigate possible 
formation scenarios, we isolated a sample of 11 pseudo- bulges from our previous HST surveys 
with WFPC2 and NICMOS and obtained for this sample additional HST/ACS imaging data. 
The relation betwee n bulge prope r ties a nd galaxy properties for these systems has been 
discussed already by ICaroUo et al. These authors found that the bulge properties 

correlate very well with those of the host galaxy, suggesting that evolutionary processes of 
the disk may be responsible for the formation of the pseudo-bulges. In this paper we focus 
rather on the properties of the population of star clusters within the pseudo-bulges. 



We have two main goals. The first is to characterize the star-cluster systems of these 
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late-type pseudo-bulge hosts and to compare them wit h the growing body of literature fo- 
cused on spirals with classical bulges (for a review see iBrodie &: Straderl |2006| ) . Under the 
assumpti on that globular clu sters form in association with enhanced bursts of star forma- 



tion (e.g. iBeasley et al.ll2002l ). the abundance of globular clusters is expected to depend on 



whether a bulge is formed "classically" , that is with a single burst of star formation, or more 
slowly, thro ugh extended star formation a s sociated with secular p rocesses in the quiescent 
galaxy disk (IKormendy fc KennicuttI |2004J : iBrodie fc Straden 120061 ). In this latter scenario, 
a lower specific frequency would be expected. 

The second goal is to continue to in vestigate the s ir nilarit ies between pseudo-bulges 
and dwarf ellipticals, first highhghted by IStiavelli et al.l (120011 ). A spiral galaxy can be 
stripped of it s disk during close interactions with neighbors, especially effective in a group 
environment (IMoore et al.lll996l ). This "galaxy harassment" scenario suggests evolution of 
late-type spirals with a pseudo-bulge into dwarf ellipticals. The stellar populations of low- 
mass galaxie s in clusters are indeed consistent with some fraction having evolved in this 
manner (e.g. IConselice. Gallagher fc Wysell2003l ). Under this scenario, there should also be 
a correlation between the properties of the star-cluster systems of these two galaxy classes. 
Here we explore this connection by measuring the specific frequencies of globular clusters 
per bulge light in our sample and comparing it with those of dwarf ellipticals. 

This paper is organized as follows. We briefly introduce the sample and discuss the data 
reduction in Section [2J Source detection and photometry are described in Section El Star 
clusters are identified in Section H] based on a color and luminosity selection. Systematic 
uncertainties are discussed in Section [51 Specific frequencies for the star cluster systems are 
derived in Section [61 while Section [71 briefly illustrates the properties of nuclear star clusters. 
Our interpretation and conclusions are discussed in Section [HI 



2. Sample and Data Reduction 



We use images of 11 late- type spiral galaxies obtained during the HST-ACS pseudo- 
bulge survey GO-9395 (PI: C. M. CaroUo). The same systems were also observ ed previously 
with WFPC2 and NICMOS JCaroUo et al.lll997l . ll998l : ICarollo fc Stiavellilll998h . All galaxies 
were selected on the basis of: (1) angular diameter larger than 1'; (2) regular morphological 
type; (3) redshift less than 2500 km s~^, to guarantee a high angular resolution in physical 
size; and (4) an inclination angle, estimated from the apparent axial ratio, smaller than 75°, 
to avoid strong obscuration of the nucle us by the disk. The objects of this study had a lso be 
shown to possess an exponential bulge (jCaroUo &: Stiavellilll998l : ICaroUo et al.l 120071 ). The 
galaxy coordinates, distance (assuming Hq = 73 km s~^ Mpc~^ for the Hubble constant). 
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absolute magnitude, morphological classification and galactic extinction are listed in Table 
m For each galaxy, we have ACS/WFC images in the F435W and F814W filters, WFPC2 
images in the F606W filter and NICMOS images in the F160W filter. 

The exposures in the F606W (V^) filter were acquired between 1996 and 1997 by WFPC2 
with the galaxy nucleus centered on the PC camera (the only WFPC2 chip considered in this 
study as its resolution is comparable with ACS/WFC), which has a scale of 0.046" pixel"^ 
and a field of view of about 35"x 35". The observations were carried out in fine lock with 
a nominal gain of 15 electrons/DN and a total exposure time per field of 600 s, split in 
two to allow removal of cosmic rays. For these images the raw data were processed with 
the standard WFPC2 pipeline (CALWP2) in order to use the most recent reference frames 
for fiat-fielding, bias and dark current subtractions; the cosmic-rays were removed using the 
STSDAS task CRREJ and the remaining hot pixels were removed by interpolation. 
Finally, sky subtraction was performed by determining the sky values from the WF chips, 
in areas farthest from the nucleus. 

The images in the F435W (B) and F814W (/) filters were acquired between 2002 and 
2003. The observations were carried out with the ACS/WFC channel, whose pixel scale is 
0. 05" pixel" ^ for a tota l field of view of 3.4'x 3.4'. Details of these observations can be found 



m 



CaroUo et al.l (120071 ). The standard ACS pipehne (CALACS) was used to perform the 
basic data reduction (fiat-fielding, bias and dark subtractions and removal of the overscan 
regions), then the IRAF STSDAS tasks ACSREJ and DRIZZLE were used, respectively, to 
achieve cosmic-ray rejection and to correct for geometric distortions. Finally, any remaining 
hot pixels were removed by interpolation and sky subtraction was carried out for all galaxies. 

The images in the F160W (H) filter were acquired between 1997 and 1998 by NICMOS 
(Camera 2). The pixel scale is 0.075" pixel"^ for a total field of view of about 19.2"x 19.2". 
These observations were carried out in snapshot mode and split into multiple exposures (to 
allow better cosmic-ray rejection) for a total exposure time per field of 384 s for ESO 498G5, 
ESQ 499G37, NGC 1345, NGC 1483, NGC 3259 and NGC2758 and a total exposure time per 
field of 256 s for NGC 406, NGC 2082, NGC 3455, NGC4980 and NGC 6384. These images 
were reduced through the standard pipeline software (CALNICA). The various NICMOS 
anomalie s (e.g. the ped e stal a nomaly) were corrected on a case-by-case basis. We direct the 



reader to ICaroUo et al.l (120021 ) for further details. 



^IRAF is distributed by NO AO, which is operated by AURA Inc., under cooperative agreement with the 
National Science Foundation. 
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3. Photometry 

3.1. Source Detections & Aperture Photometry 

Our goal is to identify star-cluster candidates and to characterize their integrated pho- 
tometric properties. Therefore we aim at constructing a uniform catalog of sources, focusing 
on detections that lie in the common area of the F435W, F606W and F814W images. The 
resulting search radius for each galaxy is given in the last column of Table [TJ The NIC- 
MOS F160W images have a more limited field of view and we use them only to evaluate the 
possible impact of dust extinction (see Section [5TT!) . 

Our study is limited to star clusters close to the nucleus of the host galaxy, where 
there are significant inhomogeneities in the image background, mainly due to dust lanes 
and gradients in the integrated-light profile. Therefore, we identify star-cluster candidates 
on a 'variance-normalized' image, that allows improved detection of point-like sources in 
regions of rapidly varying background (jMiller et al. 1997 ). The 'variance-normalized' image 



is obtained by summing the F435W and F814W images to improve the signal to noise ratio 
and then by dividing the result by a smoothed version of itself, realized with a 11x11 pixel 
median filter. 

We use the IRAF task DAOFIND to search for star clusters in the images. This task, 
optimized for point-source detection, is adequate for our goal, due to the compact size 
(FWHM ^ 1 or 2 pixels) of the star clusters in our sample. Indeed, for our purposes, with 
our data, D AOFIND marginally ou tperforms extended-source detection packages such as 



SExtractor flBertin &: Arnoutslll996l ). 



The coordinates of the sources detected on the 'variance-normalized' image are then 
mapped into the local pixel coordinates of each image using the IRAF tasks XY2RD, RD2XY 
and GEOMAP. We visually inspected all sources, with particular attention paid to those with 
FWHM > 0.25" in the F435W band, corresponding to a half-light radius R = 0.15" after 
deconvolution with the PSF, if the source has a Plummer surface brightness profile (see 
Section 13.21) . The angular scale R = 0.15" corresponds to R ~ 13 pc at 17 Mpc — the 
distance of our closest galaxy — and R ~ 27 pc at 38 Mpc — the distance of the farthest 
galaxy. These more extended sources have in general very low surface brightness and often 
appear to be artifacts and are thus excluded from the subsequent analysis. This cut is 
primarily motivated by the size distribution of galactic GCs, which have radii smaller than 
about 20 pc. In addition, all the sources in our final catalog have been visually inspected to 
remove bright foreground stars identified from their diffraction spikes. 

We perform circular aperture photometry of all the detected sources using an aperture 
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radius of 0.15" for each image and estimating the background level in an annulus between 
0.5" and 1" in radius. We verify that this method yields a correct estimate of the background 
including the galaxy luminosity. We checked this by choosing the central part of the galaxies, 
where their luminosity gradient is steepest and co mparing our photometry o f the nuclear 



-f 



star clusters against the published photometry of ICaroUo et al. I (119971 . mm . which took 
into account a detailed model of the galaxy light profile. In addition, the aperture correction 
and error estimates have been performed via Monte Carlo simulations on the actual images 
(see Section [X^ . therefore any residual bias is corrected to first order. 

We consider only sources that in all bands have a magnitude greater than a chosen 
threshold magnitude, corresponding to a signal-to-noise ratio (SNR) greater than 5 (see 
Table [2]). The photometric calibration is done by converting instrumental magnitudes to 
the VEGAMAG magnitude system by applying the photometric zer opoints listed i n Tabl e 



[21 All magnitudes are then corrected for Galactic extinction following ISchlegel et al.l (119981 ). 
Since the HST filters differ from the Johnson- Cousins band-passes, reddening corrections are 
computed using the effective transmission curve for each HST filter, created with the IRAF 
task CALCBAND in the SYNPHOT package (see Table E]). 

In this paper we indicate with uppercase letters (e.g. V) the apparent magnitude in the 
HST filters, while we use symbols like My for the corresponding absolute magnitude. Abso- 
lute magnitudes in Johnson-Cousins bands are instead indicated with symbols like Mvjq- 



3.2. Aperture Correction & Uncertainty Estimates 

Total magnitudes are determined by correcting for flux outside the measurement aper- 
ture. Aperture corrections, photometric uncertainties and the completeness of our observa- 
tions are estimated by means of Monte Carlo simulations. We add artificial star clusters 
of known luminosity and size to each image in random positions within the search area de- 
fined in Section 13.11 and we measure the retrieved flux with the same aperture photometry 
procedure adopted for the real sources. The luminosity profile of the simulated star clusters 
follows a projected Plummer law: 



where b is the assigned scale length and L is the assigned total luminosity. This is then 
convolved with the Point Spr ead Function (PSF) in the appropriate filter obtained with 



R 



TinyTim (IKrist &: Hookll200ll ) before being added to the image. 



The result for a representative photometric Monte Carlo simulation is shown in Fig. [TJ 
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Note that the distribution of the output magnitude [rriout) at a fixed input magnitude (mj„) 
has a width much smaller than the difference rriout — rriin. We define the aperture correction 
as the difference between and the median of the distribution of the output magnitude 
^out- We find that the aperture correction is almost independent of m^, and depends mildly 
on the size of the artificial source, with variations of about 0.2 mag for clusters with observed 
half-light radii between 0.05" and 0.15". We therefore apply the aperture correction to the 
observed sources according to their apparent size, using linear interpolation across our grid 
of models. 

The error on the photometry is similarly estimated from the Monte Carlo simulations 
as the variance of the output magnitudes of the artificial sources. The typical uncertainties 
are plotted in Fig. [2] as a function of the magnitude for three galaxies representative of our 
sample, NGC 1483 {D = 17.1 Mpc), NGC 4980 {D = 23.9 Mpc) and NGC 2758 {D = 31.3 
Mpc). Note that the photometric error so obtained is larger than the formal statistical 
photometric error reported by the IRAF aperture photometry task. 

We neglect the effect of possible varying charge transfer efficiency (GTE) because the 
observations were carried out when the instruments were still young, thus the effect is not 
severe. Further, the background levels in the images are provided by the und erlying galaxies. 



rather than just the sky, and are high enough to make the effect minimal (IHoltzman et al. 



19951). 



The errors on the colors are also quantified using our Monte Garlo simulations. We 
place and recover artificial sources with a given input color {Mb — Mj = B — I = 1, 1.5, 2; 
Mv - Mi = V - I = 0, 0.5, 1; Mb - My = B - V = 0, 0.5, 1) on the corresponding images. 
From the recovered photometry we measure the error on the color and correlate it with the 
errors in each band, under the assumption of a linear relation. For example, in case of the 
B — I color the error-relation reads: 

'^B-i = (rl + a]-2-rBi-crB-ai (2) 

where tbi is the linear correlation term. The typical uncertainties on the colors are plotted 
in Fig. [3] as a function of the magnitude and the color itself for the three reference galaxies 
introduced above. Our analysis yields that the typical value of the correlation r is only 
mildly dependent on the color itself, with variations from r ^ 0.5 to r ^ 0.9. We assume a 
reference value of r ^ 0.7 (the average value of the sample), which turns out to describe the 
color errors satisfactorily. 

Finally, through the same Monte Garlo simulations we also determine the completeness 
of our observations. For this purpose we define as 'successfully recovered' a synthetic source 
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whose output magnitude satisfies: 

rriin < mout < rriin + rricorr + 2(7 (3) 

where mj„ is the assigned input magnitude and rricorr and cr are respectively the aperture 
correction and the one-sigma photometric error previously estimated. Note that (i) rricorr » 
a (for example, see Fig. [1]), thus a source which is recovered as brighter than its input 
luminosity has an error at several standard deviations (essentially this happens only when a 
synthetic source falls on top of another existing source), which implies failure in the recovery; 
(ii) the error distribution is highly skewed (see Fig. [T]). Thus the completeness is similar if 
we were to change the cutoff in Eq. [3] to either la or Sex. 

Since we found that the completeness does not significantly depend on the radius of 
the objects, we show in Fig. |3]only the completeness curves derived for the mean size of the 
stellar clusters (0.075") for three galaxies representative of our sample and we note that in 
all cases the incompleteness becomes severe only at > 25.5. 



Identification of star clusters 



To classify star clusters, we compare the colors of the sources in the final catalog with 
those exp ected from synthetic stella r populations. For this we use the spectral ener gy distri- 
bution of iBruzual fc CharlotI (120031 ) models based on the Padova (1994) tracks, a ISalpeter 
(119551 ) initial mass function (IMF) with masses between 0.1 and 100 Mq, a range of metal- 
licities from 0.02 to 1 Zq and ages from 1 Myr to 15 Gyr. In addition we have also included 
sel f-consistently Hydr ogen and Helium recombination lines as well as metal lines, as described 



m lOesch et al.l (120071 ). The integrated model colors for the HST filters of the observations 
have been obtained by processing the synthetic spectra through the IRAF task CALCBAND 
within the SYNPHOT package. With the CALCBAND task we also include optioii a l dust 
extinction {E{ B — V)^l mag) . For comparison we also created iBruzual fc CharlotI (120031 ) 
models using a IChabrier I (12003! ) IMF with the same mass cutoffs. As expected, we verified 
that the resulting optical colors are identical in both cases (in fact the two IMFs are different 
only for stars below 1 M0). The mass-to-light ratio is a factor ^ 1.8 smaller (with very low 
dependence on age and metallicity) when using a Chabrier IMF instead of a Salpeter IMF, 
due to the lower fraction of low-mass stars. 

The Salpeter-IMF tracks in the — / vs. B — V color-color plane are shown in Fig. O 
and compared with the observed sources. Based on these colors, most of the detected sources 
are consistent with being relatively young star clusters (age ~ 50 Myr). The data-model 
comparison for all galaxies, except NGC 6384, suggests that only a modest extinction is 
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present (but see section 15.11 for a more detailed discussion). The spread of the observed 
sources in the cofor-color plane (standard deviation and central point of the distribution 
reported in Fig. is comparable with the photometric errors on the colors (see Fig. [3]). 
Hence it is not possible to infer the star-cluster formation history from this plot. 

A better diagnostic is provided by a color (B — I) vs. magnitude (M/) diagram (see fig. [6] 
where the Mj magnitude of the model tracks is computed for a reference 1 Mq mass). We 
can estimate the total mass of each observed source based on the difference in Mj magnitude 
compared to the model with the same B — I color. The B — I color also highlights that up to 
50% of the total detected sources are younger than 8 Myr (see Table [3]) under the assumption 
of no dust (again, see section \5A] for a discussion of the possible impact of dust extinction). 
In all galaxies, the observed sources are clearly not distributed as an evolutionary sequence 
at constant mass. Instead the older the stellar population, the greater is its luminosity (and 
hence mass, because the mass-to-light rat io increases with age). This is consistent with the 



infant-mortality scenario for star clusters (IFall et al.ll2005l ). in which only a small fraction of 



star clusters (usually the most massive) survives over time. 



Every galaxy in our sample has a number of young star clusters, with ages < 8 Myr and 
masses of the order of 10'^ M0. These can be used to reconstruct the recent star formation 
history of the parent galaxies. Taken at face value, our results indicate an average recent star 
formation rate, in such clusters, of 0.1 — 0.01 Mq/jt. Linking this SFR to the actual SFR 
is however challenging, both because we need to extrapolate our results to young clusters 
below the detection threshold and because the photometric uncertainties might introduce a 
Malquist bias difficult to quantify. Conservatively we can consider our results a lower limit 
to the overall galaxy SFR. 

In this paper we focus on the properties of the oldest (age ^ 250 Myr) and most 
massive (M ^ 10^ Mq) stellar clusters; for reasons discussed further below, we identify these 
as young candidate globular clusters. Our choice in age is an operational choice to avoid 
(i) the youngest systems, where stellar evolution of the most massive stars can influence 
the dynamical evolution of the system by their high mass-loss rates and (ii) a knee in the 
luminosity color relation (see Fig. [H]) that would lead to a degeneracy in the mass estimate 
for the sources. Once a metallicity for the stellar population is assumed, the age and mass 
selection translates into a single color and luminosity selection respectively. If we assume a 
low metallicity {Z = 0.02 Zq), the nuclear star cluster age turns out to be longer than the 
Hubble time, suggesting that these galaxies are more metal enriched (at least in their central 
regions studied here). For a metallicity range from 0.2 to 1 Zq the total number of sources 
selected as globular cluster candidates varies at the 15% level at most. In fact, the tracks at 
different metallicities are close to each other for objects older than 250 Myr (see Fig. [7]). 
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Therefore, in the following we adopt solar metallicity for the clusters. The age selection 
^ 250 Myr then translates into a color selection {B — I) ^ 0.73 mag, assuming no dust. 
The same color cut corresponds instead to a younger age selection (age > 100 Myr) if 
dust is present (see section [3?T]) . Our mass selection (M ^ 10^ M©) translates into a cut- 
off luminosity that depends on the color of the source: all the sources whose magnitude 
is brighter than the corresponding 1 Mq BC03 track, shifted by 12.5 mag, are accepted. 
Our selection has been carried out with the assumption of a Salpeter IMF. The luminosity 
cutoff we use refers instead to a mass M ^ 5.6 ■ 10^ Mq if we were to assume a Chabrier 
IMF. We identify t hese older clusters as young candidate globular clusters since we calculate 



( iTrenti et al.ll2007l ) that they have more than a 50% probability of surviving tidal dissolution 
over a Hubble time, assuming they orbit within a point-like potential well. In addition, the 
age selection we apply (age > 100-250 Myr) is sufficiently long to rule out early dissolution 
by supernova feedback. 

For each galaxy, the candidate globular cluster luminosity function (GCLF), i.e. the 
number of candidate globular clusters per unit of magnitude, is plotted in Fig. [HI For all the 
galaxies, the GCLF (completeness corrected) is reasonably fitted by a Gaussian distribution, 
whose FWHM and peak magnitude are reported in Table O The residuals from the fit 
are somewhat asymmetric with an excess of faint sources. This might refiect an intrinsic 
skewness of the GCLF, which is however difficult to assert based on our data. Alternatively, 
it might be due to the presence of a small fraction of stellar contaminants with a luminosity 
distribution peaked at the faint end of the candidate globular cluster luminosity distribution 
(see Section [52]) • Note that dust extinction does not introduce skewness in the GCLF, but 
only changes its peak value unless the amount of dust present correlates with the luminosity 
of the star clusters. 



5. Systematic Uncertainties 

The results in section [H are derived under the assumption of no dust extinction present 
outside the Milky Way and of absence of contamination in our sample. Here we discuss these 
two issues. 



5.1. Dust Extinction 



All the galaxies in our sample have NICMOS F160W coverage in the central region, 
which we use to quantify the impact of dust extinction. The F160W band is ideal, since 
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in combination with the other bands we consider, it significantly increases the wavelength 
baseline. However the more limited field of view of NICMOS does not allow us to apply this 
diagnostic to the complete sample. 

For the subsample of star cluster sources within the NICMOS field of view we perform 
a minimum chi-square fit on the data in four photometric bands comparing them to single 
stellar population models with a variable amount of dust extinction. We have three free 
parameters: the total stellar mass of the source, its age and its dust content. The fit is 
performed assuming different extinction laws within SYNPHOT: Galactic, LMC and SMC. 
As expected, the best fitting models with dust tend to be younger than their counterparts 
with no extinction. The estimation of the total mass of the sources is instead not strongly 
modified by allowing this additional degree of freedom. 

We summarize the results of these fits in Table HI where we give the number of sources 
older than 100 Myr and more massive than lO^M©, for different assumed extinction laws. 
These results are compared, in the same table, to those obtained assuming an age ^ 250 
Myr and no dust, and it may be seen that the estimates are within a factor two of each 
other. The difference between an age of 100 and 250 Myr is not very significant within the 
context of infant star-cluster mortality: in fact, most mortality happens within the first 50 
to 100 Myr and is connected to mass loss indu ced by the rapid stellar evolution of the most 



massive stars in the cluster (jParmentierl |2009| ) 



To take into account the effects of dust on the general sample, where NICMOS coverage 
is absent, we introduce a galaxy-dependent statistical correction fext-Uas to the total number 
of sources we define "massive" and "old" star clusters. fext-Uas (reported in the last column 
of Table [5]) is computed within the NICMOS field of view as the fractional difference between 
the average number of sources identified in the presence of dust vs. the number identified 
under the no-dust assumption. This correction for dust extinction is then applied when we 
use the number of old star clusters obtained from the no-dust scenario to derive their specific 
frequencies (see Section |6]). 

The correction fext-bias is conservative. In fact, based on the likelihood ratio test, the 
statistical significance of dust-extinction is above the 90% confidence level only for about 
one third of the sources with NICMOS coverage. This means that for some of the sources 
the apparent detection of non-zero extinction might simply be a consequence of allowing an 
additional degree of freedom in the fit. Finally, the NICMOS field of view i s located at the 



cente r of the galaxy, where the amount of extinction is maximal (e.g. see iHolwerda et al. 



20051 ). The star clusters outside the NICMOS field of view are thus expected to be less 



influenced by dust. 
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5.2. Sample Contamination 



A second systematic source of uncertainty that needs to be evaluated is the presence of 
stellar contaminants in our sample. These could either be stars in the Milky Way too faint 
to be identified by the presence of diffraction spikes or super-giants in the host galaxy. 

To quantify the impact of Milky Way interlopers, in each galaxy of our sample we 
selected a region of the ACS field as far away from the galaxy as possible, with an area com- 
parable to that used to search for star clusters. On these external regions we selected sources 
using the B and I bands according to the selection criteria we apply to old star clusters. The 
surface density of sources identified in the outskirts of the field of view is an upper limit to 
the number of Galactic contaminants. For all the 11 galaxies we obtained surface densities 
from 1 to 5 % (average 2.5%) of that of old star clusters selected in our main search region, 
therefore we can rest assured that Galactic stars are a negligible source of contamination 
compared to the larger uncertainties related to the treatment of dust extinction. 

To quan ti fy th e impact of bright stars in the host galaxy, we consider the models of 
Marigo et al.l (120081 ). that provide luminosities for giant-branch stars. Even for very low 



metallicities {Z = 10 '^Zq), the maximum luminosity of giant stars is My > —7; for higher 



metallicities the peak luminosity decreases (IMarigo et al.ll2008l ). From our Fig. [H] it is imme- 
diate to see that at most a few percent of the sources we select as star clusters are fainter 
than My = —7. Brighter contaminants are possible if they are young, extremely massive 
stars. T hese rare hyper giants can reach My ~ — 10 and exhibit a wide range of colors 



[e.g. see Ide Jagerl Il998l for a review). However, they are very short lived (a few million 
years) and thus not likely to be found outside the HII regions where they were born. Ideally 
we would need images in a narrow filter centered around the H-a emission line to identify 
these sources. As this is not available for our galaxies, we estimate their impact on the 
sample based on the measured star formation rate. Assuming a Salpeter IMF, we expect 
fror n a few to ~ 10 stars with mass M > 60 Mq (bright enough to reach My < —7, e.g. 
see IStothers fc Chin! Il999l ) to enter in our selection for each galaxy. Combining these two 



sources of stellar contaminants, that is Galactic and extragalactic stars, we assume a fraction 
10% of sample-contamination (indicated as fcont in Table [5]). While this might carry some 
uncertainty, its impact is certainly secondary compared to the one induced by extinction, 
except for ESQ 498G5 and NGC 6384. 
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6. Specific Frequency 



For each galaxy, we calculate the specific frequency of the candidate globular clusters 
in our field of view, i.e. the older, more massive cluster population normalized to the host 
galaxy luminosity (see Table [6]). This is defined as follows: 

-.f.OA(Mv,„ ,+15) 



s 



Ngal 



(4) 



where Myj^^^^ is the absolute magnitude of the galaxy within the search area in the Johnson- 
Cousins V-band and Ntot (reported in Tabled]) is the total number of candidate GCs, obtained 
by integrating the completeness-corrected GCLF and applying a correction fcont = 0.1 for 
the estimated fraction of stellar contaminants and a galaxy-dependent correction fext-Uas 
to take into account dust extinction (see section Wf2\\ . The luminosity of the host galaxy 
within the star-cluster-search area in the HST-F606W filter is calculated through aperture 
photometry. The background level is estimated using the WF chips of WFPC2, which cover 
a larger region of the sky compared to the star-cluster-search area, limited within the PC 
chip. The F606W magnitude {V) is then converted to the Johnson- Cousins Mvj^ absolute 
magnitude using: 

Mvjc = V + 0.287 ■{¥-!) + dm, (5) 



where dm is the distance modulus, listed in Tableland the color l^ — / is from lCaroUo et al. 



(120071 ). except for NGC 1345 (for this galaxy the color is calculated from the photometry of 
the galaxy search area in both images). The coefficient 0.287 is obtained by fittin g with a 
linear law the relation between Myj^ and V — I for the iBruzual fc CharlotI (120031 ) models 
of old star clusters (magnitudes comput ed with the IRAF task CALCBAND). Eq. [5] agrees 
very well with the calibration given by iHoltzman et al.l (119951 ) in his Table 10 (maximum 
difference of 0.05 mag. in the range — 1 < — / < 2), but has two advan tages. First it is 



calibrated exactly on the color range we are interested in (see the caveats in iHoltzman et al. 



19951 on the limited range of validity for his conversion) and second our formula is calibrated 
on the class of sources in w hich we are interested (star clusters, i.e. adopting an IMF) rather 
than individual stars as in IHoltzman et al.l (Il995l ). In computing the specific frequencies we 
have excluded NGC 1483 (the closest galaxy in our sample) because its bulge is very extended 
spatially and it gives the dominant contribution to the light within the star-clusters search 



area. 



The resulting values of S^gai are listed in Table [6l where uncertainties on Sn take into 
account the errors on the candidate GC count (Ntot), including those derived from dust 
extinction (see section I^TTI) . but not those in the distance modulus as they are negligible 
compared to the other errors. 



The definition of specific frequency adopted here is a generalization of the standard 
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Harris &: van den Berghl (jl98ll ) definition, which is based on the total luminosity of the 
galaxy and on the total number of GCs. We note though that we do not expect a large 
difference in the two definitions, because the PC camera of WFPC2 contains a large fraction 
of the total light of the galaxies. Indeed, the absolute visual magnitude Myj^gai we measure 
in our sample differs by less than about 1 magnitude compared to the total Myj^ luminosity 
of the galaxy (estimated from the absolute blue magnitude reported in Table [H assuming 
B-V = 0.5). 



We also compute the specific frequency relative to the bulge luminosity (SNtui)^ defined 



as: 



S 



Nbul 



(6) 



where again Myj^^^^ is expressed in the Johnson- Cou sins system, obtained as above from the 
bulge V^-magnitude published by lCaroUo et al.l (120071 ). The resulting frequencies of candidate 
globular clusters per bulge light are also shown in Table [61 

These frequencies for ou r sample of galax ies can be compared with published values 



for nucleated dwarf galax ies (IMiller et al 



1998) 



as well as for early and late type galaxies 
with classical bulges (see iBrodie fc Straderll2006l ). The Sjqgai we measure in our late type 



galaxies is qualitatively consistent with that of other spiral galaxies, but there is a moderate 
tendency toward higher Sn- This finding seems to be at odd with predictions from semi- 
analytics models of globular clusters for mation which link th e specific frequency with the 
overall star formation rate of the galaxy (IBeasley et al.ll2002l ): based on their scenario we 
would have expected a lower for our sample compared to that of spiral galaxies with 
classical bulges. In fact, we expec t a lower star format ion rate in pseudo-bulges if they form 
from the inner disk. However the iBeasley et al.l (|2002| ) study was aimed at reproducing the 
properties of GCs in elliptical galaxies, therefore the comparison is only indirect, based on 
the assumption that the bulge formation process is similar to that of small ellipticals. In 
addition, it did not include evolution of the GCs system of the simulated galaxies, but rather 
its properties were fixed and fine tuned to match the observed data at the time of birth. 

Another possibility to reconcile with these theoretical expectations is that of disruption 
of our candidate globular clusters as they age. Most of the 'old' star clusters in our sample 
appear in fact to be relatively young compared to Galactic globular clusters (ages from a few 
hundred Myr to one Gyr), so this suggests that tidal interaction s with the parent galaxies 
will reduce the number of star clusters as they age. For example, iGnedin &: Ostrikerl (119971 ) 
estimate that for the Milky Way, the Galactic globular cluster system has an half life of 
the order of the Hubble time. This is consistent with detailed N-b ody simulations of the 
dynamics of star clusters in the presence of a tidal field (e.g. see iTrenti et al.l l2007l and 



references therein). Under a scenario of significant adult mortality for star clusters, the 
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specific frequency for our sample will become lower than that of spirals with a classical 
bulge, as expected on the basis of theoretical modeling of star formation (SF): in the case 
of pseudo-bulges, with a likely extended but low-efficiency star formation history, fewer 
globular clusters are fo rmed than in an equivalent b urst of SF, such as that considered to 



create classical bulges (IKormendy fc Kennicuttll2004l ). 



Mean specific frequencies for nucleated dwarf galaxies ((•S'Tv^^ar/) = 7.5) — whose GCs 
are similar to candidate GCs in our sample — lie instead between our {S^gai) = 2.0 and 
{SNbui) = 30. In addition, the trend of frequencies of candidate globular clusters per bulge 
light versus luminosity is very similar to that of dwarf galaxies (see Fig. [9]). This is suggestive 
of a stripping scenario in which dwarf galaxies originate from the bulges of late -type spirals. 



One possibility is that galaxy-galaxy fiy-by encounters (galaxy harassment: iMoore et al. 



19961 ) may strip away stars in the disk, leaving preferentially behind bulge members and the 
central star clusters, both more protected because they sit deeper in the potential well of 
the system. In fact, passive evolution of the bulges (aging them from 1 to 10 Gyr, to model 
dwarf galaxies that have older stellar populations than pseudo-bulges) and accounting for 
tidal dissolution and stripping of star clusters (50 % tidal dissolution and 50% to 75% 
stripping) shift the candidate GC specific frequency vs. luminosity relation to match that 
valid for dwarf ellipticals (see Fig. [9]). 



7. Nuclear Star Clusters 



m 



A detailed study of the photome tric properties of nuclear star clusters has been presented 



CaroUo et al.l (Il997l . 1 19981 . |2002| ). here we briefiy discuss their inferred ages and dust 



content. A summary of the photometry for the central sources is reported in Tabled while 
the resulting mass, age and dust content from our maximum likelihood fit are in Tabled The 
nuclear star clusters in the sample are typically the brightest clusters of their host galaxies, 
with an inferred mass up to ~ 10'' Mq. The mass of the central star cluster broadly correlates 
with the bulge luminosity (and therefore mass), as observed in a wide sample of gala xies, 



both photometrically (e.g. IWehner fc Harris 



dynamically from spectrosc opy (e.g. 



2006 



Geha et al. 



Cote et al 



2002 



2006 



Walcher et al 



Rossa et al.l l2006l ) and 



2006). The correlation 



is consistent with Fig. 1 in Wehner fc Harris (j2006 ). Like in the Wehner fc Harris (2006) 



work, we identify a few clusters (NGC 406, NGC 1345 and NGC 2758) with a bright bulge 
10*'^ Lq), but with a relatively small nuclear star cluster mass (~ 2 ■ 10^ Mp, 



From the nuclear star cluster colors we infer a typical age > 1 Gyr in most of the sample, 
even after accounting for a variable dust extinction as described in section I5.1[ Two galaxies 
(NGC 1483, NGC 2758) present a clear evidence that their nuclear star cluster is composed 
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of a young stellar population (age 5 — 50 Myr) associated with dust {E{B — V) > 0.2). 
The nuclear star cluster of NGC 3259 has instead colors that are poorly fitted by both 
models with and without dust extinction: no reliable constrain on the stellar age can be 
obtained. The younger stellar ages observed in the nuclear clusters of NGC 1483 and NGC 
2758 might be due to a recent burst of star formation at the center of these galaxies that has 
rejuvenated them. This is not surpris i ng as rejuvenation is observed in a large number of 



nuclear star clusters ( Rossa et al.ll2006l : IWalcher et al. 



20061) and suggests a in situ formation 



scenario for the cl usters, as p r opose d by iMilosavljevid (120041 ) and inferred from observations 
of local spirals by lSeth et al.l (120061 ). However the relative old ages (> 1 Gyr) of most of the 
sources are also consistent with a formation scenario driven by mergers of star clusters that 



reach the center of th e galaxy by dynamical friction (ITremaine et al.lll975l : iLotz et al.l 12001 



Andersen et al.ll2008l ). 



8. Discussion and Conclusion 



In this paper we studied the properties of (relatively) old star clusters in a sampl e of 11 
late- type spiral galaxies selected for the presence of a pseudo-bulge (ICaroUo et al.lll997l ). Star 
clusters were classified based on their color and luminosity through comparison of population 
synthesis models with HST photometry in three bands (ACS F435W, F814W and WFPC2 
F606W). By means of artificial source detections we estimated a completeness of ~ 50% 
down to V = 25.5. The clusters in our sample present a wide range of ages and masses, 
from young blue clusters with ages of a few tens of Myrs to an older, red population (age 
> 100 — 250 Myr). We focus on these older, red clusters and identify them as young globular 
cluster candidates. 

All these galaxies have massive nuclear star clusters with masses in the 10^ — 10^ Mq 
range. These sources are typically the brightest star clusters in their host galaxy. They have 
a relatively old stellar age (> 1 Gyr) except for two cases where a younger and dusty stellar 
population is inferred. Overall, the properties of our sample appear consistent with both 
proposed formation scenarios for nuclear star clust ers, namely merg ing of stellar clusters 
driven to the galax y center by dynamical friction (ILotz et al.l l200ll ) or in situ formation 
jMilosavlievicll2004h . 



The presence of a young star clusters allows us to set a lower limit to the star formation 
rate in the galaxies, which turns out to be 0.1 — 0.01 M0/yr. This continuous, low rate of star 
formation is consistent with the formation scenario for pseudo- bulges, postulated to arise out 
of secular processes in the disk (IKormendy &: KennicuttI 120041 ). compared to a violent burst 



of SF needed to form a classical bulge. 
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Based on the iKormendy fc KennicuttI (120041 ) discussion of the formation of pseudo- 
bulges, we would expect them to have a deficit of massive star clusters compared to spirals 
with classical bulges, but this is not what we find. On the contrar y, the specific frequencies 



{number of star clusters normalized to the galaxy luminosity — [Harris fc van den Bergh 
19811 ) for the old population is consistent, within our uncertainties, to published data for 
other spirals. This appears to be a solid result, especially since we considered only central 
star clusters (the WFCP2 high resolution detector limits our area of search to the central 
35" X 35"), normalizing the specific frequency to the galaxy luminosity within the field of 
view: in general star cluster systern s are more spatially extended than their host galaxy 
(e.g. see iDjorgovski fc MeylanI Il994l : [Jordan et al.ll2009l ). thus our specific frequencies are 



probably lower limits to the global specific frequency. 

When the specific frequency is computed with respect to the bulge luminosity we get 
even higher Sn- Interestingly the specific frequency vs. bulge magnitude SNbuii^Vjcbui) 
trend is similar to the one observed in dwarf ellipticals. Pseudo-bulges have photometric 
and kinematic properties very similar to dwarf ellipticals, thus it is suggestive that some 
dwarf ellipticals might be the result of evolution of spiral galaxies with pseudo-bulges: the 
ga laxy disk might in fact be stripped during galaxy-galaxy interactions (galaxy harassment 



Moore et al.lll996l ). In this scenario star clusters are also stripped away, but still a sizable 



number might survive compared to the more spatially extended disk. At the same time the 
pseudo-bulge survives almost untouched by stripping, protected as it sits at the center of the 
galaxy potential well. While this scenario is overall appealing, detailed numerical simulations 
are needed for a proper validation. These will be presented in a follow-up paper. 
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Fig. 1. — Distribution of retrieved magnitudes (Vout) obtained through a Monte Carlo sim- 
ulation with Vin = 23 mag (red dotted line) for the NGC 2758 galaxy in the l^-band. The 
resulting median value of the distribution is also shown (blue dashed line). Sources with 
retrieved magnitude brighter than the 2a cutoff (red dot-dashed line) and fainter than Vtn 
(red dotted line) are defined as successfully recovered. 
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Fig. 2. — Photometric uncertainties for star cluster-like sources in three galaxies represen- 
tative of our sample (NGC 2758 solid line, NGC 4980 dotted line and NGC 1483 dashed 
line). These uncertainties have been estimated via Monte Carlo simulations by placing and 
recovering artificial sources in the original images. 
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Fig. 3. — Photometric uncertainties in the color measurements relative to input magnitudes 
for representative input colors and galaxies in the sample, estimated from Monte Carlo 
simulations as in Fig. O 
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Fig. 4. — Completeness curves for our survey, determined via Monte Carlo simulations by 
placing and recovering artificial star-cluster sources with a Plummer surface brightness profile 
and half-light radius of 0.075", which corresponds to the mean size of the stellar clusters in 
our sample. 




Fig. 5. — Color-color diag rams (B — V vs. V — I) for star clusters in each galaxy, showing 
a comparison of our data to iBruzual CharlotI ([2003) models. Black dots are the data with the 
"old" (age ^ 250 Myr) and massive (mass ^ 10^ M©) clusters highlighted in bold. The large red dot 
corresponds to the nuclear star cluster of each galaxy. The models include emission lines and have 
solar metallicity, single stellar population, a lSalpeteii (jl955l ) IMF and no reddening [E{B — V) = 0, 
blue line] or a 1 mag of reddening [E{B — V) = 1, green line]. Along the theoretical tracks we 
mark the age of the synthetic population in Myr. The star cluster average colors with their la 
dispersions are shown in red. 
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Fig. 6. — As in F ig. [5] but for the col o r- mag nitude diagram {B — I vs. Mj) for each galaxy 
compared to the iBruzual fc CharlotI (120031 ) models with a Salpeter IMF. Here the tracks 
for the models are shown for a total model mass of 1 Mq. Therefore the difference in Mj 
magnitude between an observed source and the model with a corresponding B — I color gives 
the total mass of the observed star cluster. In addition, we show, as dashed lines, tracks 
without emission lines. The upper envelope of the data points is a selection effect given by 
the magnitude threshold of the survey. 
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Fig. 7. — Color-magnitude diagram {B — I vs. M/) tracks for iBruzual fc CharlotI (120031 ) 
models for different metallicities (solid black line and dots: Z©; red dotted line and squares: 
0.4 Zq] blue dashed line and triangles: 0.2 Zq). The plot highlights that there is relatively 
little metallicity dependence. Ages in Myr are superimposed on the tracks. 
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Fig. 8. — Solid lines: Candidate Globular Cluster Luminosity Functions in the V-band, 
completeness-corrected, with the best-fitting Gaussian superposed. Dotted lines show the 
GCLF without completeness correction. No extinction or contamination correction has been 
applied. The black dot corresponds to the nuclear star cluster of each galaxy. The histograms 
) are normahzed so that the total number of observed sources N = Jj^^ ^^^^^c/My. 
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Fig. 9. — Candidate globular cluster specific frequencies in our sample (normalized to the 
galaxy luminosity: blue triangles; normalized to the bulge luminosity: red circles). For 
comparison the gl obular cluster specific frequencies for the sample of dwarf ellipticals of 
Miller et all fll998h are plotted as black squares. The solid red (specific frequency normalized 
to bulge luminosity) and black (specific frequency for dwarf ellipticals) lines represent a best 
fitting linear relation in the {Log^S^)', Myj^) space. The orange region delimited by the red 
dotted lines is a passive evolution of the solid red hne, assuming that the bulge population 
ages from 1 to 10 Gyr (corresponding to a dimming of ~ 2.2 mag in My) and a reduction 
in the number of star clusters due to tidal dissolution (50%) and stripping due to galaxy 
harassment (50% to 75%). 
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Table 1. Basic parameters for the sample galaxies. 



Name 


a(J2000) 


(5(J2000) 


Btot 




dm 




Type 


E(B-V) 


^search 




(h m s) 


(deg ' ") 


(mag) 


(Mpc) 


(mag) 


(mag) 




(mag) 


(kpc) 


ESO 498G5 


09 24 41.10 


-25 05 33.0 


13.96 


37.7 


32.88 


-18.92 


SXS4P/SBbc 


0.107 


3.254 


ESO 499G37 


10 03 42.10 


-27 01 39.0 


13.15 


17.7 


31.24 


-18.09 


SXS7*/SBc 


0.075 


1.529 


NGC 406 


01 07 24.10 


-69 52 35.0 


13.10 


19.8 


31.48 


-19.38 


SAS5*/Sc 


0.024 


1.708 


NGC 1345 


03 29 31.70 


-17 46 42.0 


14.02 


19.1 


31.41 


-17.39 


SBS5P/SBa 


0.038 


1.654 


NGC 1483 


03 52 47.60 


-47 28 39.0 


13.00 


15.1 


30.90 


-17.90 


SBS4/Sb-Sc 


0.007 


1.308 


NGC 2082 


05 41 51.20 


-64 18 04.0 


12.62 


17.1 


31.17 


-18.55 


SBR3/SBb 


0.058 


1.481 


NGC 2758 


09 05 30.80 


-19 02 38.0 


13.46 


31.3 


32.48 


-18.02 


PSB.4P?/Sbc 


0.126 


2.707 


NGC 3259 


10 32 34.68 


-1-65 02 26.8 


12.97 


24.8 


31.97 


-19.00 


SXT4*/SBbc 


0.015 


2.140 


NGC 3455 


10 54 31.20 


-f 17 17 02.8 


12.87 


19.8 


31.48 


-18.61 


PSXT3/Sb 


0.033 


1.708 


NGC 4980 


13 09 10.20 


-28 38 28.0 


13.19 


23.9 


31.89 


-18.70 


SXTlP?/SBa 


0.072 


2.063 


NGC 6384 


17 32 24.42 


-1-07 03 36.8 


11.14 


22.4 


31.75 


-20.61 


SXR4/SBbc 


0.123 


1.934 



Note. — Right Ascen sion (a), declination (5) and total apparent blue magnitude {Btot) are from RC3 catalog 
Vaucouleurs et al.ll99ll) . Distance (Z3) and distance modulus {dm) are from the online NASA/IPAC Extragalactic Database 
(NED). The absolute magnitude Mg^^^ is obtained from Btot and dm. Mor phological clas s ificati ons are from the RC3 {left) 
and from the UGC (right) catalogs. The color excess E(B-V) is taken from lschleeel et al.l l ll998h . The last column lists the 
radius of the common ACS/WFC-WFPC2/PC field of view. 

''using Ho = 73 km Mpc~^. 



Table 2. Details of the photometric calibration. 



Filter 


Zeropoint 


Threshold 


Reddening 




(mag) 


(mag) 


(mag) 


F435W (ACS/WFC) 


25.779=^ 


26.5 


1.319 


F606W (WFPC2/PC) 


22.084*= 


26 


0.908 


F814W (ACS/WFC) 


25.501=" 


26 


0.586 



Note. — Zeropoints used for calibrating the magnitudes, threshold 
magnitudes corresponding to a SNR ~ 5 and reddening coefficients 
expressed as Ahst filter/^Vjc ■ 

''Derived from IPavlovskv et al ] l|2006h . 

''Derived from iHoltzman et al.l \l99^ . 
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Table 3. Number of observed Star Clusters. 



Name 


All 


Age < 8 Myr 


Age > 250 Myr 


ESO 498G5 


477 


225 


62 


ESO 499G37 


434 


230 


53 


NGC 406 


1058 


376 


133 


NGC 1345 


490 


194 


41 


NGC 1483 


1208 


551 


59 


NGC 2082 


1401 


576 


96 


NGC 2758 


777 


326 


141 


NGC 3259 


805 


409 


70 


NGC 3455 


1317 


623 


67 


NGC 4980 


1070 


591 


133 


NGC 6384 


153 


25 


71 



Note. — For each galaxy (first column) the total number 
of sources identified as star clusters is given in the second 
column, while the number of young (age ^ 8 Myr) and old 
(age ^ 250 Myr) star clusters is in the third and fourth column 
respectively. 



Table 4. Star Cluster counts with extinction. 



Name 


(iV ^ 100 Myr)G„(__B^t 


{N > 100 Myr)LMC-Ext 


(N ^ 100 Myr)sMC-Ext 


(N > 250 Myr)NO-Dust 


ESO 498G5 


7 


3 


6 


5 


ESO 499G37 


2 


1 


1 


5 


NGC 406 


17 


14 


17 


23 


NGC 1345 


10 


8 


12 


18 


NGC 1483 


1 


1 


1 


10 


NGC 2082 


8 


3 


7 


18 


NGC 2758 


15 


11 


13 


20 


NGC 3259 


7 


5 


6 


14 


NGC 3455 


3 


3 


3 


7 


NGC 4980 


7 


5 


8 


19 


NGC 6384 


6 


5 


5 


6 



Note. — Number of star clusters more massive than lO'"* Mq and with age 100 Myr identified in the central region of the 
galaxies within the NICMOS F160W field of view. The number of sources has been obtained using a least chi-squared fit of the four 
band photometry (F435W, F606W, F814W & F160W) allowing for a variable amount of dust extinction with different extinction 
laws (Galactic; second column, LMC third column and SMC fourth column). The last column reports the number of sources within 
the same field of view that are older than 250 Myr and more massive than 10® Mq when the fit is forced to have no dust extinction. 
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Table 5. Candidate Globular Cluster Luminosity Function fitted as a Gaussian. 



Name 




FWHM 




fcont 


fext — bias 


Ntot 


ESO 498G5 


-9.32±0.09 


1.46±0, 


,14 


67±10 


0. 


,1 


< 0.05 


.7+10 

3'-10 


ESQ 499G37 


-7.69±0.10 


1.46±0, 


.15 


53±9 


0. 


,1 


0.70 


14+10 

^^-13 


NGC 406 


-7.86±0.09 


2.15±0, 


,15 


147±15 


0. 


,1 


0.30 


00+20 
^^-31 


NGC 1345 


-8.56±0.14 


1.80±0, 


,26 


33±6 


0. 


,1 


0.45 


16+'' 


NGC 1483 


-6.87±0.13 


2.01±0, 


,22 


61±9 


0. 


,1 


0.90 


5+4 


NGC 2082 


-8.16±0.12 


2.13±0, 


,20 


84±10 


0. 


,1 


0.65 




NGC 2758 


-8.75±0.08 


2.08±0, 


,13 


171±17 


0. 


,1 


0.35 




NGC 3259 


-8.45±0.08 


1.49±0, 


,14 


72±10 


0. 


,1 


0.55 


00+13 
^y_i9 


NGC 3455 


-7.91±0.12 


2.00±0, 


,19 


69±10 


0. 


,1 


0.55 


97+12 

^'-18 


NGC 4980 


-8.12±0.08 


1.92±0, 


,14 


139±14 


0. 


,1 


0.65 




NGC 6384 


-8.34±0.11 


1.93±0, 


,19 


76±11 


0. 


,1 


0.10 





Note. — Completeness corrected GCLF gaussian fit for the galaxies in our sample 
(first column). Fitted peak Mv-magnitude and FWHM are in the second and third 
columns respectively. The fourth column gives the total number of candidate globular 
clusters as inferred from the fit {Nf^^ ). The fifth column contains the estimated fraction 
of stellar contaminants in our sample (see section 15.2^ . The sixth column gives the 
fractional bias in the number of old star clusters induced by dust extinction of their host 
galaxy (see section|5!TJ. The last column reports our fiducial number Ntot of candidate 
globular clusters obtained, defined as Ntot = Nfit ■ (1 - fcont) ■ (1 - fext-bias)- 
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Table 6. Specific frequencies of candidate GCs. 



Name 


-^search 


^ JC gal 






^Nbul 




(kpc) 


(mag) 




(mag) 










1 4+0-2 


-Li .UD 


S fi+1-5 
^-0-1.6 


ESO 499G37 


1.529 


-17.08 


2 1+1-5 


-12.79 


110+'''5 
^I'-'-lOl 


NGC 406 


1.708 


-18.01 


^••^-1.8 


-16.65 


20t^ 


NGC 1345 


1.654 


-17.68 


1 4+0-6 


-17.57 




NGC 2082 


1.481 


-18.20 


1 4+0-7 
'-■^-1.0 


-15.65 




NGC 2758 


2.707 


-18.87 




-17.20 


13« 


NGC 3259 


2.140 


-18.55 


1 -,-1-0.5 
-'-•-'--0.7 


-14.97 


or, -1-13 
■^^-20 


NGC 3455 


1.708 


-18.01 


1 7+0.8 
-1.1 


-13.66 




NGC 4980 


2.063 


-18.29 


■^■'--1.5 


-17.39 




NGC 6384 


1.934 


-19.56 


9+°-2 


-19.44 





Note. — Candidate globular cluster specific frequencies with respect to 
the galaxy luminosity (fourth column) and to the bulge luminosity (sixth 
column). The relevant luminosities used here are given in the third and 
fifth columns and are expressed in the Johnson-Cousins system. The galaxy 
luminosity has been obtained by considering the light within a circular region 
of radius Rsearch (second column), corresp onding to the search area for the 
star clusters. The bulge luminosity is from lCaroUo et al. 1 1I2OO7I) . 



Table 7. Photometry of Nuclear Star Clusters. 



Name 


a(J2000) 


<5(J2000) 


B 




V 




/ 




H 






Rb 






(h m s) 


(deg ' ") 


(mag) 


(mag) 




(mag) 


(mag) 




( 


arcsec 





ESO 498G5 


09 24 40.68 


-25 05 31.7 


20.93±0.04 


19, 


,72±0, 


,02 


18.93±0.17 


17, 


,23±0, 


,33 


0, 


,04±0, 


,02 


ESO 499G37 


10 03 41.69 


-27 01 38.1 


20.2±0.1 


19, 


,32±0, 


,05 


18.8±0.25 


18, 


,35±0, 


,35 


0, 


,25±0, 


,05 


NGC 406 


01 07 24.48 


-69 52 30.2 


22.29±0.06 


21, 


,22±0, 


,04 


20.47±0.14 


19, 


,15±0, 


,28 


0, 


,04±0, 


,02 


NGC 1345 


03 29 31.66 


-17 46 42.6 


21.7±0.1 


20, 


.91±0, 


,07 


20.36±0.19 


18, 


,88±0, 


,31 


0, 


,05±0, 


,02 


NGC 1483 


03 52 47.65 


-47 28 37.3 


22.02±0.04 


21, 


.54±0, 


,06 


21.02±0.07 


19, 


,76±0, 


,19 


0, 


,04±0, 


,02 


NGC 2082 


05 41 50.92 


-64 18 02.6 


21.79±0.07 


20, 


,39±0, 


,02 


19.48±0.18 


17, 


,76±0, 


,23 


0, 


,05±0, 


,02 


NGC 2758 


09 05 31.17 


-19 02 33.1 


20.97±0.06 


20, 


.29±0, 


,02 


19.82±0.15 


18, 


,62±0, 


,26 


0, 


,08±0, 


,05 


NGC 3259 


10 32 34.67 


+65 02 25.9 


20.54±0.08 


19, 


.06±0, 


,05 


18.96±0.24 


17, 


,05±0, 


,33 


0, 


,04±0, 


,02 


NGC 3455 


10 54 31.12 


+17 17 04.1 


20.78±0.05 


19, 


,16±0, 


,01 


18.86±0.13 


17, 


,35±0, 


,30 


0, 


,05±0, 


,02 


NGC 4980 


13 09 10.18 


-28 38 33.6 


21.1±0.03 


20, 


.07±0, 


,01 


19.44±0.08 


18, 


,23±0, 


,11 


0, 


,05±0, 


,02 


NGC 6384 


17 32 24.27 


+07 03 36.1 


20.6±0.08 


18, 


.97±0, 


,07 


17.87±0.13 


15, 


,98±0, 


,27 


0, 


,06±0, 


,03 



Note. — Summary of the photometric properties of nuclear star clusters for our sample of galaxies, listed in the 
first column. The second and third columns report the coordinates of the central star cluster, followed by apparent 
magnitudes [B, V, I, H). The last column is an estimate of the half-light radius in the B band, after deconvolution 
with the PSF. 
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Table 8. Properties of Nuclear Star Clusters. 



Name 




Mass 


Age 


E{B - V) 


ESO 498G5 


5 - 


12 ■ 


lO^'Mg 


5-13 Gyr 


0.06-0.2 


ESQ 499G37 


1.7 - 


- 3.9 


■ IO^^Mq 


1 - 3 Gyr 


0-0.15 


NGC 406 


1.8 - 


- 3.9 


■ IO^^Mq 


2 - 5 Gyr 


0-0.2 


NGC 1345 


1.1 - 


- 1.7 


■ IO^Mq 


1-2 Gyr 


0-0.25 


NGC 1483 


0.2 


- 1 ■ 


IO^^Mq 


8-50 Myr 


0.2-0.25 


NGC 2082 


8.6 - 


- 15 


• IO^'Mq 


5-12 Gyr 


0.1-0.17 


NGC 2758 


6 - 


31 ■ 


lO^M© 


5-50 Myr 


0.5-0.7 


NGC 3259 


2.5 - 


- 21 


• 10«Mq 


5-3000 Myr 


0-1.17 


NGC 3455 


1.1 - 


- 2.8 


■ IO^'Mq 


2-10 Gyr 


0-0.24 


NGC 4980 


8.1 - 


- 26 


■ 10«Mq 


3-13 Gyr 


0-0.03 


NGC 6384 


8 - 


20 ■ 


1O^M0 


2-13 Gyr 


0.2-0.8 



Note. — Properties of the nuclear star clusters derived from the 
fit of photometric properties (see Tab. [7]l using single-stellar pop- 
ulation models. The range in Mass (second column), Age (third 
column) and Extinction (fourth column) is based on the range of 
acceptable 1 a fits. 



